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On '. ABSTRACT 

D ' We report the detection of the A1560 and A1657 ground-state C I absorption 

^ '. features in the Zabs = 1.6748 system toward the QSO PKS 1756+237. We find 

^^ no associated C I* lines with a resulting 3 a excitation temperature upper-limit 

of Tex < 8.54 (+0.65, -0.56) K, which is consistent with the predicted CMBR 

K^ , temperature of Tcmbr = 7.291 K. Because the redshifted CMBR populates the J = 1 

ly^ ' level and leaves little room for additional local excitation through either collisions or 

^SJ ■ UV fields than existing Copernicus observations of similar column density sightlines 

^^ ■ in the Milky Way. We also detect several Ni II lines and the weak Fe II A1611 line. 

Q^ . From the Ni/Fe column density ratio we find evidence for dust at a dust-to-metals 

r~| . ratio similar to that seen toward warm Galactic disk clouds. Based on these findings 

♦^1 and supported by our Lya spectrum we propose to reclassify this system as a damped 

O ' Lya absorber. 

^- 

C^ 

Subject headings: cosmic microwave background — quasars: absorption lines 
k> I quasars: individual (PKS 1756+237) — ISM: abundances — dust 



1. Introduction 

Absorption spectra of QSOs sample the distribution and physical state of intervening gaseous 
material in the early universe. Damped Lya systems (log A^(HI)^ 20.3 cm~^) are particularly 
important since they are believed to originate within young galaxies or galaxy building blocks. The 
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simple ionization state and high column densities within these systems enable the determination 
of accurate abundances for many species. Metallicities are commonly based on Zn II and S II 
measurements, while Fe-peak elements, eg. Cr II, Fe II, and Ni II, are used to infer depletion 
estimates of dust grain formation. Abundances of a-process elements, eg. Si II, Mn II, and Ti II, 
may reflect the importance of high-mass stars to the enrichment of the interstellar medium (ISM) 
within the absorbers. The excitation of C I fine-structure levels depends upon the gas pressure. 
Also, redshifted Cosmic Microwave Background Radiation (CMBR) photons excite C I more 
efficiently than the present-day CMBR due to the predicted temperature dependence on redshift, 
T = TcmbrII + z). Therefore, measurements of C I excitation in damped Lya systems can 
yield sensitive upper-limits on the temperature of the CMBR in the past, testing the fundmental 
prediction of simple big bang cosmologies that the universe cools as it expands. Neutral species 
like C I are associated in the Milky Way ISM primarily with dense, molecule-rich lines of sight. 
Only three high-redshift C I absorbers are known ( Pe et al. 1997 , Songaila et al. 1994, Meyer et 



al. 1986, Blades et al. 1985, Blades et al. 1982), consistent with the apparent lack of molecules in 
many damped Lya systems. In this Letter we report the unexpected detection of C I absorption 
sX. z = 1.6748 toward the QSO PKS 1756+237. 



2. Observations and Analysis 

PKS 1756+237 (my = 18.0, Zem = 1-721) was observed May 30 and 31, 1997 using the W. M. 
Keck Observatory I 10-m High Resolution Echelle Spectrometer (HIRES, |Vogt et al. 1994 ) for two 



hours. The signal-to-noise ratio varies from --^ 5 at the blue end to ^ 25 at longer wavelengths. 
Two slightly offset grating tilt angles yielded complete spectral coverage from 3190 - 5055 A 
with a resolution of 35,000 (FWHM = 8.5 km s~^). Full details of the spectra extraction and 
reduction procedures, along with a complete analysis of all the observed spectral features, will be 
presented in our companion paper ( Roth et al. 1999| ). Among the surprises found in the z = 1.6748 



absorption system are several Ni II lines, Fe II A1611, and C I AA1657,1560. None of these were 
expected to be present because of the low H I column density in this absorber, log A^(HI) ~ 
19.5 dTurnshek et al. 1979 ). Figure |l| shows these lines plotted versus heliocentric corrected rest 



vacuum velocity along with C II* A1336 and Si II A1527. We have examined the entire spectrum 
for other absorption systems which may confuse individual features in the z = 1.6748 absorber. 
The only possibility is a system z = 1.7337, probably associated with the QSO, which could be 
contaminating our C I A1560 line profile with Si II A1527 absorption. However, we discount it as 
insignificant because this absorber has a high ionization state and the corresponding C II A1335 
line, expected to be much stronger than Si II A1527, is quite weak (~ 20 mA). 

In this Letter we largely restrict our analysis to the weak features in Figure || which trace the 
highest column density neutral cloud harboring C I. Equivalent widths iyV\) are measured by direct 
integration across the line profiles. Statistical W\ errors are estimated from the signal-to-noise 
ratio in the adjacent continuum ( Jenkins 1973| ). To account for continuum placement uncertainties 



we inflate the errors cr(Wx) to equal 150% of the formal statistical values. The resulting rest 
equivalent width measurements and column densities resulting from a linear curve-of-growth 
(-^linear) are given in Table |^. The good A^unear agreement among different absorption lines 
from the same species suggests that most of these weak features are unsaturated. Therefore our 
column density measurements will be largely insensitive to the adopted model line profile. Our 
three-component Voigt profile absorption model is given in Table |2[ We have convolved this 
model with a Gaussian instrumental 6- value of 5.1 km s~^ to produce the superimposed profiles 
presented in Figure |l|. The total model column density in all three velocity components is given 
in the last column of Table ||. Note that the C I absorption, which presumably arises from the 
densest material, is present only in the central narrow velocity component. 



3. Abundances and Depletions 

Fe II, Ni II, and Si II are dominant ionization stages in neutral H I material and thus measure 
gas-phase abundances. Fe and Ni are both Fe-peak elements and deplete readily onto dust grains 
in the Milky Way ISM. Si is formed via the a-process and is much less susceptible to dust, but its 
abundance relative to Fe and Ni may depend on star formation history. For these reasons, the 



interpretation of relative abundance ratios is complicated (eg. Vladilo 1998, Pettini et al. 1997a 



Lu et al. 1996| ). The cosmic abundances for Si, Fe and Ni have relative values of 1.00 : 0.91 : 0.050 



(Anders &: Grevesse 1989) while our observed relative abundances (Table ffl) are 1.00 : 1.16 : 0.031. 



The observed 2 a enhancement of Fe over Si is probably not significant since the Si II A1808 feature 
is 75% deep and therefore possibly saturated. However, the observed deficit of Ni with respect 
to Fe is more likely to be real. The Fe abundance arises from the very weak Fe II A1611 feature 
which is unaffected by saturation. Therefore the Fe II column density uncertainty is simply the 
14% equivalent width measurement error. We estimate the Ni II column density 1 a uncertainty 
to be lb 3.6% from a weighted average of the six linear curve-of-growth 1 a values. This leads 
to an observed Ni/Fe abundance ratio of 0.027 ± 0.004, significantly below the cosmic Ni/Fe 
value of 0.055. Since both Fe and Ni are Fe-peak elements, this observed deviation in the relative 
gas-phase abundances from the cosmic value is likely due to depletion onto grains. Diffuse clouds 
in the Milky Way ISM come in roughly three categories with different physical conditions and dust 
depletion properties. In order of increasing temperature and decreasing density, these are cool disk 
clouds, warm disk clouds and halo clouds with gas-phase Ni/Fe values of 0.012 - 0.035, 0.028 - 
0.035, and 0.026 - 0.046, respectively ( ^embach &: Savage 1996|) . Our measured Ni/Fe abundance 
ratio is completely consistent with galactic dust depletion, perhaps suggestive of warm disk gas. 
Unfortunately, our data do not extend to long enough wavelengths to cover Zn II and Cr II. These 
lines should be easily detectable, and the relative abundance ratios Ni/Zn, Fe/Zn or Cr/Zn show 
a more pronounced variation with ISM environment because Zn is observed to deplete very little, 
if at all, onto dust grains (eg. Pettini et al. 1997b , Roth fc Blades 1995 , Sembach et al. 199^ ). 



4. Carbon I Excitation 

We report the very surprising detection of C I ground-state absorption in the z = 1.6748 gas 
cloud (Figure |2|, Table |I|). Since the lines are somewhat saturated, we have applied a 6 = 2.5 km/s 
Gaussian curve-of-growth to our equivalent width measurement errors to estimate the uncertainty 
in our model column density value (Table H), yielding iV(CI) = 1.41 (+0.30, -0.25) x 10^^ 
cm~^. We do not detect any C I* excited-state lines. While in the Galactic ISM the dominant 
C I excitation mechanism is neutral particle collisions, at high redshift the CMBR will make a 
non-negligible contribution to the population of the J = 1 level. The amount of pumping is given 
by the Boltzmann equation: 

.,_,,,. ^{,„(^)}-' „ 

where fo^i is the frequency of the energy separation between the ground and first-excited states, 



(/ii'o^i/^b) = 23.595 K (|Nussbaumer fc Rusca 19791) , and the present-day CMBR temperature 



2cMBR = 2.726 K ([Mather et al. 1994[) . Replacing the number density ratio nj=i/nj=o by the 
observable quantity Nj=i/Nj=o in Equation |l[ we predict an excited-state column density of 
iV(Cr) = 1.67 (+0.35, -0.30) x 10^^ ^,^^-2 q^ Figure | we superimpose the predicted C I* 
absorption profiles. As can be seen, the detection limit of our data is just above the C I* line 
strengths produced by the redshifted T = 7.291 K CMBR. 

To improve our C I* detection limit, we have shifted the observed data points to a common 
velocity scale for each of the four excited-state features (for the C I* A1561 doublet we use 
an /-value weighted average line center). This average C I* spectrum is shown in the inset 
to Figure y, where we have superimposed the predicted C I* profile produced using the same 
stacking method. The rms error of the inset stacked spectral region (S/N = 45) corresponds to 
an equivalent width error over nine pixels of 0.75 mA. We inflate this aw^ by 50% for the same 
reason as before and adopt a la uncertainty in the stacked C I* equivalent width of 1.1 mA. 
The stacked C I* profile has an equivalent width of 2.1 mA, which is consistent with no local 
excitation of the C I fine-structure levels above CMBR pumping at the 2 a level. However, since 
inspection of the stacked C I* spectrum suggests that some absorption might be present, although 
at somewhat less than our 2 a level, we have elected to adopt a 3 a detection upper-limit for this 
stacked C I* feature of Wx{Cl*) < 3.3 mA, which implies a 3o" CI* column density upper-limit 
of A^(CI*) < 2.67 X 10^^ cm~^. This leads to a 3 a upper-limit on the C I excitation temperature 
in the z = 1.6748 absorber toward PKS 1756+237 of Tex < 8.54 (+0.65, -0.56) K. The stacked 
profile for this 3 a detection limit is shown by the dashed line in the inset to Figure 0. 

In addition to measuring the CMBR temperature, our C I fine-structure population 
measurement allows us to determine the physical conditions within the absorbing gas cloud. For 



collisional excitation by hydrogen atoms, we employ the rate coefficients reported by Launay &: 



Reueff 1977. For UV pumping, we solve the detailed balance population equations for the 200 



strongest C I UV lines (de Boer & Morton 1974, Jenkins et al. 198^ , Morton 1991, Zsargo et al, 



1997) and measure the UV field in units of the WJl average interstellar radiation field ( [Morton 



p^ [Table 9], |de Boer et al. 1973| , [Witt fc Johnson 1973| ). In the absence of CMBR pumping and 



collisional excitation, we find a 3a upper-limit on the UV field within the z = 1.6748 absorber 
of < 14.3 (+3.2, —2.6) times the Galactic UV radiation field. Alternatively, if collisions alone 
populate the J = 1 level and the gas kinetic temperature is 10 K, typical of cool disk clouds, the 
hydrogen atom density is {3 a upper-limit) n^^ < 950 (+1079, —369) cm~^. The C I collisional 
excitation rate is very sensitive at low temperatures. For example, at a slightly higher temperature 
(15 K) the density limits are reduced significantly {nj^^ < 191 (+64, —44) cm~^). At kinetic 
temperatures similar to less dense clouds or halo clouds (100 or 1000 K) the upper- limits are 
further decreased to n^°^ < 25.1 (+5.8, -4.6) and n^oooK < io.6(+2.4, -1.9) cm'^. 

These limits become more interesting when one includes the considerable contribution to C 
I* excitation by the high redshift CMBR. As a result, observations of C I fine-structure lines in 
QSO absorbers can yield far more sensitive probes of the physical conditions within distant gas 
clouds than similar observations within the Milky Way can provide. For example, if we assume 
the CMBR temperature at z = 1.6748 is given by the predicted big bang value of 7.291 K, we 
find that our UV field and hydrogen atom density upper-limits are substantially lower. With 
CMBR pumping, the ambient UV field must be less than (3o") 5.6 (+3.3, —2.6) Galactic fields, 
njj°^ < 371 (+652, -224) cm-^, njf^ < 75 (+54, -37) cm'^, nJj°OK < 9.8 (+5.8, -4.6) cm'^, and 
^lOOOK < 4.1 (+2.4, -1.9) cm-^ 



In a survey of 27 Galactic stars using the Copernicus satellite, Jenkins et al. 1985 did not 
report an excitation temperature lower than 40 K and excited-state lines were detected in 23 lines 
of sight. However, excitation temperatures were not determinable for nearly 25% of their sample, 
and all four of the sightlines lacking C I* absorption had A^(CI) values that were comparable to or 
below that found in the z = 1.6748 absorber. Because of the redshifted CMBR, our limits on the 
physical conditions within the z = 1.6748 absorber are 2-3 times lower than the values in these 
four Milky Way sightlines. It is perhaps too early to tell whether high-redshift C I absorbers are 
similar to Milky Way low column density sightlines, or if there is a deficit of high column density 
C I absorbers toward QSOs. The discovery of more high-redshift C I QSO systems will eventually 
answer the latter question, while higher quality C I spectra of the Milky Way ISM are required to 
address the former. 



5. A Damped Lya Absorber? 

The presence of weak, singly- ionized metal lines and C I argue strongly that the z = 1.6748 
absorber toward PKS 1756+237 has been erroneously classified as a moderately low column 
density system and we have assumed it is in fact damped. Our spectral data extend far enough 
into the blue to cover the H I absorption feature, and we have attempted to verify the validity of 
this claim. The Lya line, smoothed and rebinned to 20% the resolution, is presented in Figure ^ 
with model Voigt profiles for the originally reported logA^(III) ~ 19.5 ( p?urnshek et al. 197S| ) and 



a damped Lya absorber superimposed. The data are noisy, the Lya absorption is located near 
the blue wing of the QSO Lya emission (zem = 1-721), and it extends across two echelle orders, 
all of which complicate the continuum placement. For these reasons, we do not claim a conclusive 
H I column density measurement. However, we believe there is sufficient evidence to reclassify 
this system as a damped Lya absorber with ~ 6-10 times more neutral material than previously 
reported. 
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Table 1. Rest Equivalent Widths and Column Densities 



Transition A 

(A) 



S/N 



Wx±la 

(mA) 



iVfinear ± l(T 

(cm~^) 



(cm 



model 
-2^ 



Fe II A1611 


1611.2005 


0.00102 


26 


27.9 ±3.9 


(1.19 ±0.17) X 10^^ 


1.23 


X 


IQis 


Si II A1808 


1808.0126 


0.005527 


19 


153.5 ±7.6 


(9.60 ± 0.48) X 10^^ 


1.06 


X 


10^5 


Ni II A1455 


1454.8420 


0.05954 


17 


44.6 ±5.5 


(4.00 ± 0.49) X 10^3 


3.30 


X 


10^3 


Ni II A1317 


1317.2170 


0.1458 


6 


83 ±15 


(3.71 ± 0.67) X 10^3 








Ni II A1370 


1370.1320 


0.1309 


9 


92 ±12 


(4.23 ± 0.55) X 10^3 








Ni II A1710 


1709.6000 


0.06884 


25 


58.7 ±4.7 


(3.30 ± 0.26) X 10^3 








Ni II A1742 


1741.5490 


0.1035 


24 


84.5 ±5.1 


(3.04 ±0.18) X 10^3 








Ni II A1752 


1751.9100 


0.06380 


26 


65.5 ±4.8 


(3.78 ± 0.28) X 10^3 








C I A1657 


1656.9283 


0.1405 


25 


29.2 ±2.4 


(0.855 ± 0.070) X 10^3 


1.41 


X 


10^3 


C I A1560 


1560.3092 


0.08041 


23 


22.1 ±3.1 


(1.28 ±0.18) X 10^3 




// 




C IP A1336 


1335.7077 


0.1149 


7 


92 ±12 


(5.07 ±0.66) X 10^3 


1.10 


X 


1014 



'^Oscillator strengths (/) from Morton 1991 except for Fell A1611 /-value taken from Cardelli fc 
Savage 1995| 

Column densities assuming linear curve-of-growth, A^unear = (1.1296 x l^^'^)Wx\~'^ f~^ cm~^ 

■^Voigt profile model column densities: Fe II and Ni II values result from a 3-cloud sum, C I value 
includes only a single velocity component 



Table 2. Three-Component Voigt Absorption Model 



Ion 


log(iVi)^ 


log(iV2)^ 


log(iV3) = 


CI 




13.15 




Fell 


14.58 


14.75 


14.47 


Mil 


13.15 


13.00 


12.95 


Sill 


14.58 


14.65 


14.37 


CII* 


13.29 


13.90 


13.05 



a, 



vi = — 19 km/s, 6i = 15 km/s 
^V2 = km/s, 62 = 2.5 km/s 
^V3 = +12 km/s, 63 = 10 km/s 
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Fig. 1.— Weak Ni II, Fe II, C I* and C I features along with saturated Si II A1527 in the z = 1.6748 
absorption system toward PKS 1756+237. The spectra have been normalized and the observed air 
wavelengths have been converted to vacuum heliocentric velocities, with v = km s~^ corresponding 
to z = 1.67477. Superimposed are Voigt model line profiles. 
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Fig. 2. — C I A1657 (top panel) and A1560 (bottom panel) ground-state lines at z = 1.67477 
toward PKS 1756+237. Superimposed on each line is our adopted model profile. Also plotted are 
the expected strengths for the four excited-state C I* lines with the redshifted CMBR as the only 
source of C I excitation. In the inset figure (bottom panel) the four C I* lines have been shifted 
to a common heliocentric velocity scale and averaged. Similarly, the expected C I* profile and our 
3 a upper-limit have been stacked and superimposed. 
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Fig. 3. — Lya absorption at z = 1.67477 toward PKS 1756+237. The data have been smoothed 
over 11 pixels and rebinned to a dispersion of 10 kni/s. Superimposed on the bottom panel are 
A'^(HI) = 3.16 X 10^^ and 2.0 x lO^'^ cm~^ Voigt profiles with the same 6-value and velocity found 
from the C I lines. 



